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Abstract 

The production of hadrons in Ni-Ni at the GSI is considered in a hadronic 
gas model with chemical equilibrium. Special attention is given to the abun- 
dance of strange particles which are treated using the exact conservation of 
strangeness. It is found that all the data can be described using a temperature 
T = 70 ± 10 MeV and a baryon chemical potential \i B = 720 ± 20 MeV. 
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I. INTRODUCTION 



Recent experiments at the GSI laboratory have shown a relatively large yield of K~ 
mesons |],[2| in collisions involving Ni+Ni at 1.93 GeV A. This result is surprising because 
the beam energy is below the threshold for the reaction N + N — > N + N + K + + K~ 
which needs at least 2.5 GeV in the lab. The production of K~ mesons therefore proceeds 
entirely via medium effects in the nuclear matter of the colliding nuclei. An explanation for 
this phenomenon has been suggested in several models @-|5[] wherein it is argued that the 
mass of the K~ decreases when nuclear density increases while the mass of the K + increases 
slightly. This would make it easier for K~ mesons to be produced since the threshold for 
pair production becomes considerably lower. The production of kaons is therefore a test of 
the effects of the density of the medium and should be analyzed carefully. 

It has been pointed out |||7| that the relative abundance of hadrons (excluding kaons) 
in the final state can be described by using a hadronic gas model with chemical equilibrium, 
i.e. only two parameters, the temperature T and the baryon chemical potential fiB describe 
all hadronic abundances except the kaons. 

It is the purpose of the present paper to address the abundance of kaons in detail. 
It is well-known that imposing exact strangeness conservation introduces a suppression of 
strange particles if the system has very small values of the dimensionless quantity VT 3 . 
For large values of VT 3 this suppression disappears rapidly. The size of the Ni system 
is relatively large however the temperatures involved are always below 100 MeV and thus 
these corrections have to be considered seriously. We will show in the present paper that all 
hadronic ratios, involving strange particles, can be explained with only three parameters, 
namely the temperature T « 70 MeV, the baryon chemical potential fis ~ 720 MeV and 
the radius of the system, R « 4 fm. The large variety of mesured production ratios of 
non-strange particles are explained with the same set of parameter values. 

There exists evidence that particle abundances in heavy ion collisions at higher energies 
are also very close to chemical equilibrium (for a review see e.g. ref. ||). The data from 
CERN favor a region around T « 160 — 200 MeV and a baryon chemical potential around 
fig « 180 — 350 MeV. The data from BNL favor a region having a lower temperature, 
T « 100 — 140 MeV and a larger chemical potential \xb ~ 450 — 600 MeV. The results from 
GSI are in a region which is T « 60 — 100 MeV and a baryon chemical potential \is ~ 700 
MeV. The data from GSI thus correspond to a low temperature but a considerably higher 
value of the baryon chemical potential \ib- Many questions are still left open. It doesn't 
explain where the K~'s are coming from or why they seem to be so "easily" produced despite 
being sub-threshold. The evidence for chemical equilibrium is very strong in our opinion. 

What about thermal equilibrium? Clearly the momentum spectra deviate strongly from 
a simple thermal distribution given by a Boltzmann factor, however, many of the effects 
cancel out when considering ratios of hadrons as we do here. Effects due to transverse flow 
or excluded volume cancel out in ratios (under circumstances described in ||). In this paper 
we would like to investigate the production of different particle species as observed at the 
GSI in Ni-Ni collisions. 
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II. FORMALISM 



To impose strict strangeness conservation one projects the standard (grand canonical) 
partition function, Z(T, Xb, Xs, Xq), onto the state with strangeness S 

Z s = ±- 'd<f> e- iS *Z(T,X B ,X s ,X Q ), (2.1) 

Z7T JO 

where the fugacity factor Xs has been replaced by 

X s = e*. (2.2) 

Since the incoming state has S = we will consider only this value in the rest of this 
analysis. For simplicity we present below the formalism used in Boltzmann approximation. 
It is more complicated to use quantum 'statistics but it can be done (see e.g. [p!QH T2|. The 
partition function is then composed of a sum over the different particle species. Those that 
have strangeness zero are untouched by the above operation, those that have strangeness +1 
acquire a factor e l< ^ while those with strangeness -1 acquire a factor e~ % ^. If only strangeness 
0, ±1 particles are present we have 



J S=0 



i- ^ d<t> exp {/V s=0 + M s =ie i4, + fV s =-ie-^} , (2.3) 



where J\fs=o is the number of particles having zero strangeness, i.e. mainly pions and nu- 
cleons. jVs=±i is defined as the sum over all particles having strangeness ±1 omitting the 
fugacity factor due to strangeness, i.e., the contribution of A particles to J\fs=-i would be 
given by 

Afs=-i = VJ 0r 3 e- E ^ T+ ^ T . (2.4) 

We have checked that the presence of cascade particles does not modify our final results 
for the analysis of the GSI data. Of course for higher temperatures cascade particles and 
Q's should be included. The strange and anti-strange particle numbers are not equal since 
there is a net baryon number in the system, e.g. there is a difference in the numbers of A's 
and A's due to the presence of the baryon chemical potential, 

1 r 2n 

Z s=0 = 7T d( t> ex P Ws=o + (jVs=i + jVs=-i) cos + % (j\f s =i - Ns=-i) sin 0} . (2.5) 

Z7T JO 

Despite its appearance, the above expression has no imaginary part. Exploiting the sym- 
metry properties of the integrand, this can be rewritten as 



1 f 2lT 

7j s=o 



r2ir 

/ dcj) exp {J\f s=0 + (JV s =i + Ks=~i) cos 0} cos [(J\f s=1 - J\f s =~i) sin 0] . (2.6) 
Jo 



2tt Jo 

The partition function can thus be written in a compact form as 



- = Z -^- r dcj) cos{xsin0}e ycos *, (2.7) 

27T JO 
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where 



x = X s =i - Afs- 



(2.8) 



and 



y = Af s =i + Ns=-i 



(2.9) 



and Z is the part of the partition function which contains only non-strange particles. 
As an example, the number of kaons is given by 



N K = gV 



d 3 p 
(2tt)s 



-E k /T 



\ J* 7 " d<p cos {<p + x sin <p} e y cos ' 



- Jq 77 d(f) cos {x sin 0} e y 



COS <f) 



(2.10) 



It is of interest to investigate the small volume limit of the above expressions since this is 
where the effects due to the canonical formalism are most obvious. A typical term has the 
following form 



p2n 






/ d 


f> exp 











which in the small volume limit will lead to 



^5=0 



1 + 



2tt 



2tt 



gv 



d 3 p 
(2?if 



-E/T+i4 



d 3 p 



-E/T+iq 



d 3 p 



-E/T-iq 



+ 



(2-11) 



(2.12) 



Note that each time there will be two terms since otherwise the integration over leads to 
a vanishing contribution. The number of will thus be given by 



N K + w V 



(2^)3 



<i 3 p 



+J v7 



d 3 p 
(2^)3 



-E A /T+fi B /T 



(2-13) 



This shows clearly that the number of particles increase quadratically with the volume. For 
a large system the dependence on the volume becomes linear. The additional suppression 
of strangeness can be seen clearly in the small volume limit above. In order to balance 
strangeness, the production of a kaon has to be accompanied by either an anti-kaon or by a 
strange baryon. This is explicitly present in equation |2.13| . 

The results are not changed substantially if one includes also particles having strangeness 
±2 and ±3 using the analytical methods of fll3,14|. 



III. NUMERICAL ANALYSIS 
A. Resonance Width 



As long as the baryon chemical potential is below ~ 800 MeV the Boltzmann approxima- 
tion is adequate. For higher values one has to include effects due to quantum statistics. Due 
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to the low temperatures involved the width of resonances has to be taken into account. For 
example, an appreciable number of pions is coming from the decay of A resonances below 
the mass of the A but still within the width of the A. For very small values of the width T, 
the Breit-Wigner resonance shape can be replaced by a 5 function: 



lim —-. — 

T^O 7T (s 



mT 



\2 _|_ m 2p2 



8(s — m 2 ). 



(3.1) 



For large widths this is no longer possible and one has to keep the Breit-Wigner resonance 
shape. For example the integration of the Boltzmann factor has to be replaced according to 



d 3 p exp 



y/p' 2 + m 2 
T 



d 3 p J dsexp 



\/p 2 + s 



1 



mT 



7r (s — m 2 ) 2 + m 2 T 2 



(3.2) 



The mass a/s is integrated from m — 2T to m + 2T. In some cases, e.g. A(1232), A(1440), 
the lower limit goes far below the threshold limit, e.g. for the A resonance -y/s = + m^. 
In such cases the lower limit is chosen to be the threshold value + m n . We have taken 
care of the fact that the normalization must be adjusted accordingly since the integral over 
the Breit-Wigner factor should still give unity. 



B. Isospin Asymmetry. 

In an isospin asymmetric system there are initially four parameters, namely the tem- 
perature T, fugacities A# and Xq, and the volume V. The charge chemical potential can 
be eliminated by considering the ratio of baryon- and charge content of the system. The 
baryon- and charge densities, ns and rig, and thus the corresponding chemical potentials, 
are related by the condition 

n B (T, \,R) = 2 n Q (T, A, R), (3.3) 

where (^) measures the isospin asymmetry in the system. For the Ni-Ni system considered 
here, one has 



^ Am - 1.04 (3.4) 



V /Ni-Ni ^ Ni 

and so there is a 4% deviation from the isospin symmetric case 



C. Results 

The hadronic gas model contains all particles listed by the Particle Data Group ]I5] . Our 
main result is shown in figure 1. Each experimentally measured hadronic ratio corresponds 
to a band in the (T,/j,b) plane. The width of the band corresponds to the error bar which 
has been reported. It is to be noted that several of these bands are almost parallel to the 
temperature axis. All these bands are fairly insensitive to the chosen hadronic volume. The 
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main exception is the K + /n + ratio. This ratio is highly sensitive to the value of the radius 
of the hadronic gas. In figure 1 we show this for a radius of 4 fm. The dependence can be 
seen more quantitatively in figure 2 where one can see the K + /n + ratio corresponding to a 
radius of 3.7 fm and also for 4 fm. Also shown is the large volume limit which is denoted by 
TD (thermodynamic) limit, corresponding to infinite volume. The temperature necessary to 
reproduce this ratio rapidly goes down with increasing values of the radius. This is because 
for a small volume there is an intrinsic suppression factor which makes it necessary to go to 
higher temperature in order to reproduce the same ratio. The region of overlap has a value 
for the baryon chemical potential of 

Hb = 720 ± 20 MeV. (3.5) 

It is not so easy to fix the temperature because many of the lines are almost parallel to each 
other, but a temperature interval given by 

T = 70 ± 10 MeV (3.6) 

will give a good fit to the hadronic ratios. 



D. The Ratio A/K 

A requirement of the exact strangeness conservation in an ensemble, where the hadrons 
i with \Si\ < 1 are included, leads to condition 

K + + K° + (K*) + A + £ + (¥*) = K' + K^+ (K*) + A + £ + (Y*), (3.7) 

where hadron symbols stand for their densities at freeze-out. Now the baryon chemical 
potential is quite convincingly fixed to the range 700 MeV < (1b < 800 MeV, and the 
temperature is found to be around T ~ 70 MeV, so the anti-strange baryons can be neglected 
due to condition Y/Y oc exp(2/i B /T) ~ 10 s . Further, the heavier resonances (denoted by 
superscript *) are not supposed to play any significant role in energies considered, so 

A + S K- . . 

1 " T^^r (3-8) 



K+ + K° K+ + K ' 

At the first stage, the effect of isospin asymmetry in the system is neglected, and we are led 
to the result (E° and A are not distinguished) 

A K~ . , 

1 - — = 0.96 ±0.02, (3.9) 



K+ K+ 

where the experimental result K + / K~ = 26 ± 9 is used. 

When the charge chemical potential is taken into account, equation (|3.8| ) leads to ex- 
pression 

A ~ (l + A Q 1 )-(l + A Q )f^ 

K+ l + cosh(lnA Q ) ' V ' J 

where Aq = exp(/XQ/T) is a charge fugacity. Using values 60 MeV < T < 80 MeV and 
700 MeV > fiB > 770 MeV we extract the corresponding values for the charge fugacity (-2 
MeV > fiQ > -6 MeV), and obtain the result A/K + = 0.99 ± 0.03. Similar procedure gives 
for neutral kaons the result A/K° = 0.94 ± 0.03. 
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IV. SUMMARY 



The data from GSI show very good agreement with a hadronic gas in chemical equilibrium 
having T « 70 MeV and ~ 720 MeV. When comparing the GSI results with those 
obtained at higher energies a clear picture emerges. There is an increase in the temperature 
and at the same time a decrease in the baryon chemical potential hb in going from GSI to 
Brookhaven (BNL) and then on to CERN. This is indicated schematically in figure ^| where 
we have made use of the recent compilations of Sollfrank |8| (see also [|l6fl ) and of |7|]. It is 
not trivial to translate the values for T and hb into values for the energy and the baryon 
densities. This is because one needs to know the particel numbers and not just the ratios as 
we have done in this paper. An estimate of these has been made recently in [[17]]. A good 
agreement with chemical equilibrium does not mean that the particle spectra should follow 
exactly a Boltzmann distribution since the momenta of particles can be severely affected 
by flow. As an example, a model with Bjorken expansion in the longitudinal direction will 
still have its particle ratios determined by Boltzmann factors even though the longitudinal 
distribution is nowhere near a Boltzmann distribution ||. 

The main deviation from chemical equilibrium is observed in strange particle abundances. 
These clearly deviate from chemical equilibrium but it seems that a single parameter mea- 
suring the deviation from chemical equilibrium is sufficient to describe most of the strange 
particles. The necessary suppression can be fully accounted for by using exact strangeness 
conservation and there is no need to introduce a strangeness suppression factor 75 as is the 
case for particle abundances measured at higher energies at CERN. 
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FIGURES 



FIG. 1. Curves in the (/j,b,T) plane corresponding to the hadronic ratios indicated. The 
interaction volume corresponds to a radius of 4 fm, and the isospin asymmetry is B/2Q = 1.04. 

FIG. 2. Curves in the (/Ub,T) plane showing the dependence on the radius of the interaction 
volume. Interaction volume corresponds to radius 4 fm. The K + /ir + and (f)/K~ ratios with the 
freeze-out radius R = 3.7, and the thermodynamic limit for the ratio K + /ir + are presented with 
discrete labels. 

FIG. 3. Location of the freeze-out temperature and baryon chemical potential for different 
energies. 
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